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The massless flow between successive minimal models of conformal field theory is 
related to a flow within the sine- Gordon model when the coefficient of the cosine potential 
is imaginary. This flow is studied, partly numerically, from three different points of view. 
First we work out the expansion close to the Kosterlitz-Thouless point, and obtain roaming 



behavior, with the central charge going up and down in between the UV and IR values of 
c = 1. Next we analytically continue the Casimir energy of the massive flow (i.e. with real 
cosine term). Finally we consider the lattice regularization provided by the O(n) model in 
which massive and massless flows correspond to high- and low-temperature phases. 

A detailed discussion of the case n = is then given using the underlying N=2 
supersymmetry, which is spontaneously broken in the low-temperature phase. The "index" 
trF(— 1) F follows from the Painleve III differential equation, and is shown to have simple 
poles in this phase. These poles are interpreted as occuring from level crossing (one- 
dimensional phase transitions for polymers). As an application, new exact results for the 
connectivity constants of polymer graphs on cylinders are obtained. These results and 
points of view are used in the following paper to discuss the appropriate exact S'-matrices 
and the resulting Casimir energies. 
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1. Introduction 

Prompted by conformal field theory, there has been enormous progress in the last 
several years in understanding two-dimensional field theories off the critical point. For both 
massive and massless models, much has been understood by using perturbed conformal 
field theory. In the numerous cases where the model is integrable, a number of non- 
perturbative techniques have been extensively developed and applied. There are many 
massive models whose properties are very well understood, but there exist some very basic 
massless ones whose properties have not been fully explored. 

A particular example is the flow between the (t + 1)^ and t^ 1 unitary models [3J, 
which has been studied perturbatively for large t 0] and for finite t using the conjectured 
thermodynamic Bethe ansatz relations 0. Moreover, this problem is related to an even 
less-understood model, which is the sine-Gordon model with an imaginary potential. This 
model is interesting in its own right, because at a particular coupling it describes the 
"dense" phase of self-avoiding polymers. These models are related because here and in 
many other cases, one can consider perturbed conformal field theories as restrictions of 
affine Toda theories, thus extending free-field representations of conformally invariant field 
theories. In the conformal case and the massive case, these restrictions are well known. 
The (Virasoro) minimal conformal field theories can obtained from a free boson || and 
the massive integrable perturbation of these theories by the $13 operator can be obtained 
from sine-Gordon ]5]|]]. In statistical-mechanical language, this is often called mapping 
the system onto a Coulomb gas M. The restriction is most easily understood in field 



theory using quantum-group structures [BT5 , and to some extent parallels lattice-model 



constructions |16|,|17|]. A very similar restriction applies also to the massive 0(n) model 

HQ 

We discuss a similar truncation for massless flows, where both the UV and IR limits 
are conformally invariant. In this case, since both UV and IR fixed points are obtained 
by restriction of a free boson, one expects that the full trajectory is the restriction of an 
intriguing flow in sine-Gordon, with central charge equal to 1 at both extremities. This 
flow actually was first mentioned in M, and related explicit field-theory calculations were 



done in ||20|| . We discussed it briefly in [^TJ where we presented numerical data with a 
"roaming" trajectory that violates the c q q conjecture. We describe in this paper and in 
the following our attempts to understand quantitatively this problem. 

The models describing the flows between the minimal models are unitary, but the 
sine-Gordon model with imaginary potential has a Lagrangian that is not real, and hence 
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is non-unitary. The non-unitarity allows a variety of interesting properties to appear, 
including level crossing, an increasing c-function, spontaneous breaking of supersymmetry, 
and singularities. While these features are part of what makes this model interesting, they 
also often make it difficult. In this paper and its sequel we will use a variety of methods 
to explore the situation. We believe that we have completed the picture for the minimal 
models; we will give the exact S'-matrices and ground-state energy. However, we will see 
that there are still a number of unanswered questions in the case of sine-Gordon. The first 
paper contains detailed studies (including numerical) of the problem, excepting attempts 
to build a scattering theory. The second paper is then fully devoted to that question. 

In section 2 we discuss the flow within sine- Gordon in an expansion close to the 
Kosterlitz-Thouless point. We recover the roaming trajectory with a pair of extrema of 
order 1/t 3 , and discuss the behavior of the stress-energy tensor and alternate definitions 
of the central charge. 

In section 3 we discuss in detail the lattice regularization of the problem provided 
by the 0(n) model on the honeycomb lattice flowing to the low temperature phase. A 
quantum-group symmetry present all along the flow is exhibited that will provide useful 
information about level crossings. A lattice derivation of the flow of independent sectors 
of the partition functions is also provided. This section may be skipped by the reader who 
is not interested in lattice questions. 

In section 4 we discuss numerically the massless flow using two approaches. We 



discuss first the analytic continuation of the massive flow ||21|| . The first extremum of the 
central charge agrees nicely with the expansion around the Kosterlitz-Thouless point, but 
breaks down afterward. We also consider lattice data, which show the existence of multiple 
extrema for small enough t, a feature not seen in the large-t expansion. 

Section 5 contains a short discussion of the boundary case n = —2, which is exactly 
solvable by other methods. The behavior is somewhat pathological, but the appearance of 
singularities provides compelling evidence for similar behavior to occur at least at nearby 
values of n. 

Section 6 contains a detailed discussion of the case n = 0. This case is physically 
very important because it describes polymers, and is also quite amenable by analytic 
techniques due to the underlying N=2 supersymmetry. We find that N=2 supersymmetry 
is spontaneously broken; the theorem that ordinarily does not allow this does not apply as 
a result of non-unitarity (despite the non-unitarity, the described model is meaningful for 
condensed-matter physics). It is shown that an infinite number of level crossings separate 
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the UV and IR fixed points. These crossings are studied analytically by calculating the 
index tr F(— 1) F |2^,^3[, where they show up as poles. The locations of these poles are 



found by solving the Painleve III differential equation in the complex plane. 



2. Thermal flow: generalities and the large-t expansion 

Consider the (t — 1)^ minimal unitary conformal field theory with central charge pi 

6 



t(t + l) 



2.1 



where t is an integer. As is well known, it can be obtained by twisting and reducing a 
Gaussian model J|] (a free scalar field 4>) which has coupling in "lattice-like" conventions 

<7 = i + i 0. 

We perturb this critical theory by the thermal operator, which is identified in the 
conformal language as $13, and has left and right dimensions 

h = h= 1 —^. (2.2) 
t + 1 y ' 

Perturbation by $13 corresponds in the Gaussian model to a perturbation by the operator 
cos 2(fi with weight, with respect to the untwisted stress energy tensor 

h = h=-!—. (2.3) 
t + 1 K J 

The coupling constants are different due to the difference in conformal weights in the 
twisted and untwisted theory. If for the minimal model 

Smin = Smin + / ^ d ^ ( 2 - 4 ) 

then 



Ssg 



f d 2 x \^(d<P) 2 + (^) 1 / 2 cos 2(/> ] . (2.5) 
J L4tt J 



The mass in the minimal model is proportional to l^l - ^ - , where 5(3 = f3 c — (3 refers to 
distance from the critical inverse "temperature" . 

For 5 f3 > the minimal model is expected to flow to a trivial massive phase. From the 
Gaussian point of view, we indeed have the standard sine-Gordon model. It is integrable, 
with a known S'-matrix, and only the massive soliton and antisoliton are in the spectrum 



for t > 2 [25 1 - From this S'-matrix and the thermodynamic Bethe ansatz (TBA), relations 
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can be written for the ground-state energy which restrict to the ones conjectured for the 



minimal models H when appropriately twisted [21 



When 5(3 < the minimal model flows |l],0 to another conformal field theory with 
central charge in the infrared 

c=l- , 6 . (2.6) 

The corresponding situation in the sine-Gordon theory is a little unusual. It leads in ( |2.5|) to 
a purely imaginary coupling, a non-unitary case that has already attracted some attention 
|7|, p0| , |26[| . The behavior can be understood qualitatively as follows. Call A e the coupling 
of the vertex operator e ie ^ in a general action with the Gaussian term plus perturbations. 
In ( |2.5| ) for instance we have A2 = A_2 = ^(SP) 1 / 2 . Approximate renormalization group 
equations can be written for the A e and the coupling constant g [|2~7| , |2lf1 . If b is the rescaling 
parameter, one finds at first order 

d a \ e = f2-^A e) (2.7) 

and for the coupling constant 



db V 2 f 

-^ = 8vr 2 A e A_ e . {2.1 
gab 



No magnetic charges arise in these equations. Observe that cos 2<p is relevant for g > 1, 
irrelevant otherwise. 

In the case 5(3 > 0, A2 = A_2 are real numbers. Then if initially g > 1 as in the case 
here, A2 and A_2 both increase along the renormalization group flow ( |2.7| ). This increase 



contributes to increasing g as well due to (|2.8|) and in turn to making cos 2<p even more 



relevant. So at this order, more and more charges ±2 are generated at large distance, and 
one expects (slightly changing the language) the "Coulomb gas" J^f^f to flow to a massive, 
plasma phase. 

Now suppose 5(3 < 0. Initially, the magnitudes of both A2 and A_2 both increase in the 
renormalization. However, since they are purely imaginary, their product is negative, so 
due to ( |2.8| ) the coupling decreases. This decrease means that the operator cos 2(f) becomes 
less relevant. One therefore expects that floating charges ±2 are present in large numbers 
at intermediate distances, but ultimately disappear in the IR. If this is so, the sine-Gordon 
model with imaginary coupling should flow back to a Gaussian model, although with a 
smaller coupling constant. This phenomena was first discussed in |7j in relation to the 
dense 0{n) model to be discussed in subsequent sections. 
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There is a simple way of understanding this heuristically [|5|,|6. 20|1 . This is to postulate 



that the action for the conformal minimal model in fl2.4| ) is in fact the Liouville action 
with a curvature term to represent the charge at infinity: 

Smin = / £ [W) 2 + iQR(f\ + A 2 e 2 ^, (2.9) 

where R is the scalar curvature. This model is reasonably well- under stood (at least with Q 
imaginary, in our conventions) and has a fixed point which is a conformal theory of central 
charge c = 1 — 6gQ 2 . Thus the conjecture is that with appropriate choices of g and Q 
(real here), this action describes a minimal model. The operator with dimensions of $1,3 
turns out to be exp(— 2i<p). Thus to perturb the model, we add A_2 / exp(— 2i0) to the 
action, and obtain the sine-Gordon model ( |2.5| ) plus a curvature term. Thus in this sense, 
the actions ( |2.5[ ) and ( |2.4| ) differ only by the curvature term. Calculations in field theory 
including the curvature term have been done for large t in [^D[ ; these indeed reproduce the 
flow of the central charge between successive minimal models. The sign of the perturbation 
is crucial: we can always redefine the coefficients of the exponential terms by shifting the 
field (f), but the product X2X-2 does not change. Thus if this product is positive we are in 
the massive phase, while if it is negative we are in the massless phase. 

2.1. Large-t expansion 

It is possible to study this a little more quantitatively in the limit of large t. Setting 

e = 2 in (|2.7Q and ( p78|) , both the fugacities and the renormalization group eigenvalue 

2 

y e = 2 — |- are small in this limit, and we can simultaneously expand in both variables. 
At dominant order these two equations allow us to compute 

^| = 4( 9 -l), (2.10) 
dg do 



which we can integrate to 

dg_ 
db 



(2.11) 



where the integration constant is determined by the UV value guv — 1 + j- This beta 
function vanishes at another value of g, indeed smaller than one as expected, gm = 1 — j 
(this value is in fact exact, but all we get here is a check of it at first order). The flow can 
be integrated to give 
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and 

i e 2(b-b )/t 

Both the running coupling constant g and the fugacities exhibit the behavior qualitatively 
described in . In particular there is a value of b = bo where the fugacities are extremal, 
so that their associated beta functions vanish. However, this is not an intermediate fixed 
point at this order, because the beta function for the coupling constant g is not zero but 
instead extremal while g passes through g = 1. 

Further information about this flow in the Gaussian model is obtained by considering 
the central charge. Recall that at first order the trace of the stress-energy tensor for the 



Gaussian model perturbed by vertex operators of charge ±e reads p8 



6 = -2ny e (A e e ie * + X^e**) . (2.14) 

From this and the relation between the derivative of the running central charge defined in 
and O we find the running central charge at distance R 

c(R) = 1 - 6n 2 y e (i^A e i^A_ e ) , (2.15) 

Put now e = 2. When 6(3 > the product of fugacities is positive, so c decreases in 
agreement with the c-theorem. On the other hand, when 5f3 < the c-theorem is not 
expected to hold due to the non-unitarity of the theory.0 In fact, one has the opposite 
behavior; since the couplings are purely imaginary, c initially increases above its UV value 
1, going against a "c e g-theorem" (this was observed in |[26j| ). Ultimately, c must return to 



1 in the IR. At the order to which we are working, c is given by 

12 e 4(6-5o)An - e 4(6-6o)A) 



! + TT ~ \„ , -73 \ (2-16) 



t 3 (I +e 4(6-b )A) 3 

and it exhibits a pair of extrema at e^ - ^! 1 = (2 ± \/3) with values 



c ± = l ± ^. (2.17) 



The R dependence disappears at this order. 



1 Here this non-unitarity manifests itself rather trivially by the negative sign of the two point 
function of 0. 
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As a result of the vanishing of the beta functions for the fugacities, the two-point 
functions (00) and (TO) vanish when g passes through g = 1 at this order. It is likely 
that higher-order contributions give them a finite but small value. Indeed, assuming that 
the coupling g evolves monotonically, its beta function never vanishes, and there cannot 
be any intermediate fixed point. 

Since c(R) is equal to one both in the UV and IR where the sine-Gordon model flows 
to a Gaussian model, for any value of R it must exhibit at least a pair of extrema. Using 
the Callan-Symanzik equation 

R M + ^r° (2 - 18) 

for coupling constants gi (here g : A) we see, without any reference to a given order of 
perturbation, that at these extrema the derivative has to vanish, which in turn means 
vanishing of ( 0(0)0 (!?)). However we can expect a priori that the peculiar values of the 
couplings gi at the extrema depends on R so there is no special intermediate point in 
the flow where (O(O)O(i?)) would vanish for every R. It is possible for such a two-point 
function to vanish for some value of R because of the non-unitarity of the theory. 

Beyond this order, the perturbative study of the flow looks quite complicated. How- 
ever, according to the arguments of 0, this model is integrable, and non-perturbative 
methods should be applicable. Recall indeed that in the massive flow, the running central 
charge c(mR) (m is the soliton mass) for the sine-Gordon model was successfully calcu- 
lated by using the TBA [jnj. (As we will explain in section 4, the running central charge 
is defined via the Casimir effect on a cylinder.) It can also be computed in the large-t 
expansion |2| and one finds at first order the same formula as (|2.15| ) with R = 1. The 
study of this function c(mR) for finite t, and the search for its corresponding TBA, are the 
subjects of this paper and the next. An immediate result that can be obtained for finite t is 
the exactness of gm = 1 — ~. This follows most easily from the 0(n) model discretization, 
and symmetries (see the discussion of eq. 3.20). We therefore see that the massless flow 
in sine-Gordon should be observed for t > 1 only. In the standard parametrization this 
corresponds to 4% < f3g G < 8tv. In particular we do not expect this massless flow to have 
a semi-classical version since the latter is obtained as (3sg ~~ * 0. 

In the large-i expansion we considered only A±2, but the renormalization group trans- 
formations generate higher charges as well. In the region where X±2 are found decreasing, 
they disappear for two reasons. On the one hand some pairs of opposite charges annihilate 
and contribute simply to renormalizing the coupling constant g. On the other hand pairs 
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of identical charges coalesce to increase first of all A±4. Of course these are irrelevant with 
respect to the UV fixed point. Their equation of evolution is made of a damping term 
2/4 A±4 and a source term proportional to So we expect A±4 to increase driven by the 
amount of ±2 charges, and subsequently to decrease when charges ±2 have disappeared. 
The picture generalizes to all higher charges. Therefore, provided the domains of evolution 
of these various charges are separated enough, one can expect more structure in c(mR) 
than just a pair of extrema, but it is difficult to make quantitative predictions. 

3. 0(n) vertex model: Symmetries of the thermal flow 

We discuss in this section various properties of the 0(n) model on the honeycomb 
lattice. This model provides a convenient regularization of our problem when it flows to 
the low-temperature "dense" phase. It is also useful in identifying symmetries of the flow 
and the occurence of level crossings. 

3.1. 0{n) vertex model and its symmetries 

It is well known that a convenient lattice regularization for the sine-Gordon model 
and its restrictions is provided by the 0(n) model on the honeycomb lattice (the links of 
a hexagonal lattice). The original partition function in the absence of a magnetic field is 

Z = Y,$ M n L . (3.1) 

where the summation is taken over all configurations of self-avoiding, mutually avoiding 
loops that can be drawn on the lattice, M is the number of lattice edges occupied by a 
piece of loop (monomer) and L is the number of loops. As it stands, ( |3.1| ) is not defined 
locally. A local description is obtained as follows M . We still associate a weight (3 to each 
monomer, but we do not give weights n to loops any more. Instead we put an orientation 
on each loop, together with a weight e ±l7re °/ 6 per left (right) turn. On a planar honeycomb 
lattice, the number of left minus the number of right turns of a closed loop is equal to ±6. 
Therefore setting 

n = 2 cos 7reo, (3-2) 

we see that summing over all the possible orientations reproduces the desired weight n per 
loop, now as a sum of local contributions. We refer to this model of oriented loops as the 
0(n) vertex model, or vertex model. It has been given other names and is well known 
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to be related in some limits to the Zamolodchikov-Fateev model and the Izergin-Korepin 



model; see [g9| . As we will discuss below, the correspondence between ( j3TT|) and the vertex 
model is a little more intricate in non-trivial topologies like the cylinder or torus. 

The difference between the non-local and local formulations is manifest if one tries to 
write a transfer matrix. In the first case, "past" information must be taken into account: 
we need to know if two monomers on a given column were connected at former times (see 
30(1 and references therein). In the second case, we simply have a vertex model with an 
.R-matrix acting in the product of two spin-one representations, and the transfer matrix 
is the appropriate tensor product of these R matrices. To see this consider propagation 



along one of the three equivalent directions in the honeycomb lattice. Following [[JT 
decompose the interactions into elementary three-leg vertices with weights 



we 



iv 7 = (3 \e, e,e, e \e \e 1 

(3.3) 



uj[, . . . , u'j = p \e \e \e \e,e,e, 



where we set e = e l7re °/ 6 . As usual we can modify these weights (by a so-called gauge 
transformation) and leave the partition function invariant. The general transformation of 
interest reads 



_ „_! c'e a'e b'e a' b' c' 
u>±, . . . ,a>7 — p , — , — — , — , — , — ,— 
a b c ce ae be 

, x a b c ce ae be 
c'e a'e b'e a' b' a 



(3.4) 



Let us set 



q = -e 6 , (3.5) 
and -A; = —q, jf^ = -<? _1 , c 2 = c' 2 = ab = a'b' . We get then 



, , , , a — 1 ,—3 3 3 —3 —3 ^3 

lox . . .,u}? — p ,e ,e,e,e ,e ,e 

,,/ _ a-\ ,3 ^-3 ^-3 ^3 ^3 -3 

io x . . . ,u 7 — p ,e,e ,e ,e,e,e . 



(3.6) 



Let us now stick two such three-leg vertices together. The rules of the model completely 
determine the state of the intermediate edge from the ones of the outer edges, so we can 
as well forget it to get a standard four-leg vertex, whose weights determine a .R-matrix. 
This matrix has two eigenvalues: with degeneracy 6, and /3 -2 — q — q~ x with degeneracy 
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three. The respective eigenvectors read 



and 



1,1) 

1,0) + 5 |0,1) 

l,-l) + g|-l,l); |1,-1)-9 _1 | 

-l,0) + g|0,-l) 

-1,-1) 

|io> - g-^o, 1> 



i,i)-^(i + <T 2 )|o,o) 



(3.7) 



(If 

|1, -1> - - 1, 1> - ^9 _1 I0, 0> 



(3.8) 



|0,-l)-qr|-l,0) 

Clearly this ^-matrix has a "partial" U q sl{2) symmetry, where the generators S only act 
on the states | ± 1) as in a spin one-half representation (acting on |0) they just vanish). 
In the following we set S z = ±1,0. The symmetry U q sl{2) therefore acts on the loop 
degrees of freedom without seeing the empty spaces in the system. By usual coproduct 
applications, this symmetry extends to the transfer matrix, up to boundary effects. Let us 
emphasize that the 0(n) vertex model has more symmetry than what we are pointing out 
here. In particular, at the solvable point of [|32| there is a U q >sl(2) symmetry acting on a 
full spin-one representation, with q' = (— e 6 ) 1 / 4 . In addition, in the low-temperature limit 
(3 — > oo there appears a U q »sl(3) symmetry [|33[ with q" = (— e 6 ) 1 / 2 . Notice the equality 



n = —q — q 



(3.9) 



We restrict in what follows to the case n G [—2, 2] for which q is a complex number of 
modulus one. 

Let us now discuss boundary conditions. As is well known, the correspondence be- 
tween the vertex and the 0(n) models breaks down in a non-trivial topology due to non- 
contractible loops. For instance, on a cylinder non-contractible loops in the space direction 
have a weight 2 in the vertex model because they have the same number of left and right 
turns. This can be repaired by changing the vertex-model boundary conditions. Such an 
operation is also called putting a charge at infinity, or a seam ||34j| . Draw on the cylinder 



a line in the time direction. Give a vertex configuration an additional weight e ±l7re every 
time an oriented edge crosses the line going upwards (downwards) (e is well defined modulo 
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two). In that case a non-contractible loop in space direction gets a weight 2cos7re which 
can be made equal to n by the choice e = eo- On a torus there are non-contractible loops 
with a more complicated topology. To reproduce their correct weight one has to consider 
a family of sectors of the 0(n) vertex model characterized by their "magnetic charge" 



m = S z (a topological defect S z it) and the "electric charge" e (a twist term) |35|,|36|]. The 
sector with — e and — S z has identical properties. Let us parametrize 

n = -2cos7T7, q = e i7r7 , 7 G [1, 2]. (3.10) 

so in particular 7 = (eo + 1) mod 2. The various sectors of the vertex model are then 
connected by quantum-group commutative diagrams [17]. This means that for k greater 



than S z = m, and k — m integer, all levels of the sector S z = k, e = rwy are contained in 
the those of the sector S z = m, e = k'j. 

3.2. The critical point 

The 0(n) model is critical when 0] 

/3 c = (2 + v / 2^r 1/2 . (3.11) 

Introduce solid-on-solid height variables <p dual to the arrows, with step value ±7r. In the 
continuum, these degrees of freedom become a free bosonic field with coupling constant, 
in the scale where topological defects are not renormalized, 

9 = 1- (3-12) 

The above parameters e, m are now interpreted as the usual electric and magnetic charges 
in a two-dimensional Coulomb gas, with associated conformal weights h em (h em ) = 3(7^= ± 
rriy/g) 2 . The partition function of the vertex model in the S z = m, e sector, and on a torus 
of size R x T reads in the continuum 

Z em -> Z em = e fR I V y ^ + ^K +2j , m . ( 3 . 13 ) 

v(ymy) ~ 

for the simple geometry we consider, y = exp(— ttT/R) = y. The central charge of this 

vertex model is equal to one for all values of n. 

The aforementioned coincidences of levels are easily observed in ( |3.13j ). For example, 

2 

the ground state of the sector S z = m, e = has weight h = h = ^t—, while weights in the 
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sector S z = 0, e = rrry are h = h = (g m + 2 -?) ; coinciding with the above for j = 0. Such 
coincidences are well known to occur in the continuum limit as following from duality in 
the Gaussian model. We stress that they in fact occur for the finite lattice model, and for 
all values of (3. Notice that j = is not necessarily the ground state of this sector: when 
m = 1 we have |2 — gm\ < g so that j = —1 is the ground state for g > 1. 

The partition function of the 0{n) model is obtained by proper summation of ( |3.13| ) 
for various sectors. An important fact is that the ground state of the 0(n) model is not 
the one of the vertex model. Choosing e = eo indeed gives the central charge 

9 

In conformal language, the 0(n) model is thus described by a twisted Gaussian model. In 
the 0(n) point of view, the value of c = 1 for the untwisted Gaussian theory is a c q q value 
associated with the presence of a negative-dimension operator. 

When q is a root of unity, the U q sl(2) symmetry implies that additional levels coincide. 
We restrict for simplicity to the case 

n = 2 cosy, t = 1,2,. ... (3.15) 

and restrict also to sectors S z = m an integer, where there are an even number of non- 
contractible loops in the space direction. In this case the commutative diagram of [17| 



extends further. Levels in the sector S z = m,e = k'j (with k > m) that are not in any 
other sector are selected by the restricted partition function 

^ab = ^ ^ ^e=kj,m+jt ^e=m~/,k+jt (3.16) 

jez 

where 

a = k-S z , b = k + S z . 



Differences like (|3.16| ) are exact for finite lattices, due to the U q sl(2) symmetry. 



Notice that the central charge of the 0(n) model reads in the parametrization ( |3.15| ) 



c=l--. -. (3.17) 

t(t + l) K ' 

The 0(n) model at integer t is not the same as the minimal model except for t = 2, 3 
(n = 0, 1) ||37|| . However, in the continuum limit the above restricted partition functions, 
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as discussed in [17] for a similar case (Potts model), are appropriate pieces of the minimal- 
model partition functions i 



Z a b + Za+l/~f,b+l/-y ~^ Zab — Xai 



(3.18) 



s=l 



for 1 < a, b < t — 1, where x are the usual characters of the corresponding irreducible 
representations of the Virasoro algebra for central charge ( |3.17| ), and the non-universal 
bulk free energy / has been discarded. The above expressions sum up to the ones worked 
to express the various minimal-model partition functions as linear combinations 



out m 



of Gaussian ones. For example, the diagonal (i.e. A-series) partition function reads 



Zaa = Z c 



t + 1 



t + 1 



(3.19) 



where Z c is the Gaussian partition function in the notation of [[JB]. The expression for 
finite lattices analogous to (|3.19| ) also makes sense due to the U q sl(2) symmetry. Since 
this symmetry holds in fact for any (3, all these expressions also make sense off the critical 
point. 



3. 3. Massless flow 

As discovered in |7j] the 0(n) model is also critical for (3 > (3 C , with a continuum limit 
that is /3-independent. In other words, this low-temperature phase is a flow to another 
fixed point, which turns out to also be a Gaussian model. This phase is often referred to 
as the "dense phase" by reference to polymers (n = 0). Of course the associated massless 
flow should be precisely the twisted version of the one discussed in section 2. In particular 
since the U q sl{2) symmetry is present all along the flow, and coincides in the UV and IR 
with the Gaussian model duality, the relations between n, g and 7 valid at (3 = f3 c should 
still hold true. This does not require g at the low-temperature fixed-point to be the same 
as in the massive phase; as we saw in section 2, we expect the coupling to get smaller. 
Thus another solution of these relations is applicable 0; we have exactly 



9 



7 



1 

1 - -. 

t 



(3.20) 



2 It is necessary to add sectors with k = I/7 to restrict the overall number of non-contractible 
loops to be even. 
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The corresponding central charge is given again by ( |3.14j ) which reads here 

C=1 -7(^T) <3 ' 21) 

After proper restriction ]T2| , |IT5| , |T5| ] , this flow gives rise to the massless flow from the minimal 
model at central charge ( |3.17|) to the one with ( |3.21|) . From the lattice point of view, as 
emphasized already, all coincidences of levels in the vertex model transfer matrix hold for 
every value of f3. The sectors corresponding to the 0(n) model or the minimal models can 
be extracted by appropriate combinations of vertex partition functions. One checks easily 
that in the dense phase 

t-2 

Z a b + 2 a+1 / 7)b+1 / 7 — > Z a b = ^^XraXrbi (3.22) 

r=l 

so we can see the flow of various sectors of the minimal models: 

t t-2 

^XasXbs SOWS tO ^XraXrb ( 3 -23) 
s=l r=l 

Such a result implies in particular that minimal models flow to a model of the same kind 
in the ADE classification [[38] . 



3. 4- Symmetries in the continuum 

We notice that the above U q sl(2) symmetry involves the same deformation parameter 
as the U q sl{2) symmetry exhibited in the continuum theory ( |2.5| ) since the formal proof for 



the existence of this symmetry |15[ does not depend on the sign of 6(3. We also recall that 
the soliton S'-matrix describing the massive flow exhibits this symmetry as well. To a large 
extent one can formally identify the soliton trajectories with the loops of the underlying 
lattice model thanks to Q HH- 

In the following we will use results for the continuum as well as the discrete theory, 
and usually refer to results in the latter case with an upper index u l". 



4. General numerical analysis: analytic continuation of massive TBA and 
lattice data 

4-1. Analytic continuation 

A natural way of getting insight into the phase structure of the 0{n) model is to 
determine its ground-state energy Eq(R), where space is a cylinder of circumference R. 
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There are two methods we will use to find this quantity. The first is to explicitly diagonalize 
the lattice hamiltonian (or transfer matrix — we shall refer to these two interchangeably) 
on a small, periodic lattice. The second way relies on the fact that when spacetime is a 
cylinder of radius R and a large length T, the partition function is Z ~ exp(— Eo(R)T). If 
we reverse the roles of space and time, this is the partition function for a large system at 
temperature 1/R. 

If the 0(n) lattice model were integrable for j3 7^ f3 l c one could use the lattice Bethe 
ansatz to diagonalize the Hamiltonian for a large number of sites and then do the ther- 
modynamics using the resulting equations for the energy levels [39]. Unfortunately, this 
does not seem possible here. However, this model is integrable in the continuum. Thus 
we can use a related but different method, called the thermodynamic Bethe ansatz (TBA) 
||4(J|| . Here one conjectures the exact S- matrix for the particles of the theory, and uses 
the relation between the S'-matrix and the phase shifts of the wave function in order to 
find the equations for the energy levels. Like in the ordinary Bethe ansatz case, one then 
uses these energy levels to minimize the free energy. The S- matrix for the sine- Gordon 
model (|2.5|) with 5f3 > has been known for a long time [p5[ , and the TBA equations 
were written in [^T| . As we discuss in the sequel, there are complications in understanding 
possible S'-matrices and the TBA for 6(3 < 0. It is thus useful to consider first analytic 
continuation of these high-temperature results to the massless phase. In this subsection we 
do this numerically; there seems to be no way of continuing the TBA equations directly. 
In |] a similar procedure essentially reproduced the flow for the ground state energy from 
the tricritical to the critical Ising model, a result which was then confirmed by using a 
"massless S- matrix". 

Our results can be written in terms of the dimensionless quantity mR, where m oc 
(5/3) ( t+1 )/ 4 . (The coefficient of proportionality can be obtained by comparing the results 
of conformal perturbation theory and the TBA. Similarly, on the lattice m oc (5(3 l )^ t+1 ^ 4 , 
where the coefficient is determined by numerical study.) We define the quantity c(mR) = 
6REo(R)/n : to take advantage of the fact that in unitary theories c(0) and c(oo) are the 
central charges of the UV and IR conformal field theories p3,ffT|. We note that off the 



critical points this is not the c-function defined in [|I|, but at lowest order in SP the two 
are the same. It is expected (and has been shown in all known unitary cases) to decrease 
monotonically as mR is increased. In non-unitary theories c(0) and c{oo) are not the 
central charges but rather c q q = c — 12h — 12h, where h and h are the lowest conformal 
dimensions in the sector p3,il]. Efforts to find a c e g-theorem have used this observation 
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as a starting point. This makes the ground-state energy a logical candidate, but as we saw 
in section 2, it does not always decrease! 

The TBA gives a set of integral equations which allow the determination of c(mR) in 
the high-temperature phase. They have not been solved directly, but it is straightforward 
to solve them numerically. There exists an expansion around mR = 

c(z) =c uv + c hu i k (z) +^f n z n , (4.1) 

n 

where we have defined z = (mi?) 4/, ' t+1 '. We also use the variable r = mR. This corre- 
sponds to the perturbative expansion of (|2.5| ) in powers of 5/3 (because of the Z2 symmetry 
cf) — > —((>, there are only even powers of (5/3) 1 / 2 in the expansion). The non-universal bulk 
term does not appear in perturbed conformal field theory and so is not necessarily an 
integer power of z. For t integer this bulk term is the same for the full sine- Gordon theory 
and for the minimal model, where it was computed in |J. This is because the models 
differ by boundary conditions only; the truncation does not affect bulk properties. More 
precisely one can check numerically (and analytically in related models like the 8-vertex 
model) that in the IR the levels of the vertex-model transfer matrix that belong to the 
minimal model are still at finite distance in the 1/R scale from its ground state, implying 
that all these levels share the same bulk term. 

To do the analytic continuation, we first solve the TBA equations for c in the massive 
phase to double-precision accuracy. We then subtract off the massive bulk term. For 
t even, this term vanishes, but for t odd it is proportional to (mR) 2 ln(mi?) [[§. After 
subtraction, we have a power series in z which can be continued to negative values. We 
do this in one of two ways, both involving Pade extrapolation. 

The first method is to find the ratio of polynomials which go though a set of z > 
data points (zi, c(zi)) determined numerically using the TBA integral equations. The ratio 
is unique after one specifies the orders of the numerator and the denominator; their sum 
is the number of data points plus one. Typically, the most accurate results are obtained 
by using 10-15 data points. We have used the routine in Numerical Recipes [|i2| |. It is 



particularly useful because it provides a convenient way of estimating the error. By using 
different subsets of the original data points, one obtains different extrapolations. The error 
is then estimated by looking at how the result changes for a given value of z. 

The second method is to find the coefficients f n in ([4.1] ) by fitting them to the 
numerically-determined c(mR) function, using a standard least-mean-squares program. 
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One then analytically (it is convenient to use Maple or Mathematica) finds the Pade ap- 
proximant. (For a polynomial of order p, the Pade approximant is the unique ratio of 
polynomials of order p/2 such that when this is expanded around z = 0, the first p + 1 
terms reproduce the original polynomial. If p is odd, then the polynomial in the numerator 
is of order (p + l)/2 with the denominator of order (p — l)/2). The Pade approximant 
should describe c for negative z for much larger values of \z\ than does the polynomial. 
This method has the disadvantage that it is it is not as easy to determine the errors in- 
volved, and it requires more computer time, because many more data points are needed 
for a good fit than for the first method. For example, in the t = 2 (n=0) case, one finds 

c(z) = 1. - .4454536z + .091380z 2 - .00903z 3 - .00018z 4 + .00015z 5 , (4.2) 

where the numerical error is in the last digit given. The resulting Pade approximant is 

, . 1 - .301756z + .04362 2 - .003142 3 , t . 

c(z) = = . 4.3) 

v ! 1 + .143698z + .0162z 2 v ! 

We have checked for several values of t that the results obtained by these two methods are 
consistent, but we have mainly used the first method. 

It is easy to find the first extremum of the c-function using this continuation. For large 
t this allows us to check formula ( |2.17| ) for this maximum, which is determined by quite a 
different route. Before making the comparison, we must subtract off the low-temperature 
bulk term. The situation here is a little more complicated than for the massive flow. For 
< n < 2 the bulk term is the same for the full sine- Gordon model and the minimal 
theory, and we can use the results of ||. This is not true for n = (t = 2) due to a 
large number of level crossings which we shall discuss further. For the sine-Gordon bulk 
term, our results are still consistent with putting t = 2 in the minimal-model formulas. 
The agreement is very good, as exemplified in the following table that gives the ratio of 
Cmax — 1 for the continuation and ( |2.17| ) 

t=2 t=3 t=4 t=5 t=6 t=7 t=8 

.556 .685 .759 .806 .840 .87 .88 

We certainly see the 1/t 3 behavior of the maximum being reproduced in the continuation, 
and standard extrapolation as a function of 1/t shows very good convergence to the value 
one. As t increases the number of coupled equations in the massive TBA increases, making 
the calculation more difficult, as well as the effect smaller to see. Other features of the 
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large-t expansion can also be checked against the TBA results, like the ratio of the r values 
at the maximum and at the intersection with c = 1. 

We have not been able to extend the continuation far enough to observe a minimum, 
nor c going back to one in the IR. We give an example of data obtained in figure 1 for 
n = V2. 

It is difficult to determine precisely the radius of convergence of the z expansion for 
c(mR), but several tests provide a value around z c « 1. The Pade extrapolation typically 
gives reliable data up to z ~ 3 for n = 0, and less as n increases. Singularities in the 
complex plane can be localized using the Pade approximants, but their position depends 
too much on the number of points used to be reliable. There was however no sign of these 
singularities getting too close to the real axis, negative or positive. (If there were a pole 
on the negative axis, the terms in ( fl.2|) should alternate in sign.) 

4-2. Lattice results 

A complementary tool for studying the massless flow is the numerical diagonalization 
of the 0(n) vertex model transfer matrix on cylinders of lattice width R l . (Since we will 
use results for the continuum as well as the discrete theory, we often refer to results in 
the latter case with an upper index "/".) This provides lattice data that, after proper 
extrapolation \(3 l — (3 C \ — > 0, R — > oo with r l = \(3 l — (3 l c \ v R l fixed, should reproduce 
the continuum limit of interest (r and r l being related by a non-universal proportionality 
constant). A difficulty here is to extract the bulk term. Several obvious strategies are 
possible, none of them giving too-satisfactory results. Their comparison however allows 
the determination of c(r l ) with a reasonable accuracy, generally better for small or large 
values of r l and for t large. A convenient number to consider is 

logA^i^m') _ logA (fi 2>(m') 

c^\r l ) = R[ 1 1 ^ . (4.4) 

where R\ and R l 2 are neighboring widths (eg R\ R l — 1). In the study of critical points such 
combination is known to converge very rapidly to the exact value c, with the bulk term 
conveniently extracted. In our situation, we are not allowed a priori to use this procedure, 
since apart from a regular term, the bulk term (actually its singular part from the lattice 
point of view, but that is all what is seen in the continuum limit) also scales like 1 / R l in 
the limit we are interested in. On the other hand, corrections to scaling for finite R l seem 
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to compensate conveniently so that (|4.4| ) is a rather reliable estimate, especially at t large, 
i.e. it coincides with estimates obtained by other, more correct means, up to error bars. As 
an example we give also in figure 1 the result obtained from (fOj) for n = y/2. After proper 
rescaling of the lattice and continuum variables, it is seen to nicely coincide with the result 
from TBA extrapolation in the region where they are reliable. For another example we 
give in the following table the maximum of the function of c^ R } (for R[ = R l , R l 2 =R l -l) 
for various sizes 

#'= 4 5 6 7 8 

Cmax= 1-024 1.017 1.016 1.015 1.015 

which seem to be converging nicely to the TBA value c max = 1.0137. In the lattice 
simulation (but not the TBA) one observes further c reaching a minima, and going back 
to the value 1 in the IR. Details of such curves can be checked against the large t results, 
like the value of c at the minima, the ratio between the values of r l at the maxima and 
minima, or at the maxima and the point c = 1, with very reasonable agreement. 

Two additional features must be mentioned. First we have observed a general tendency 
for the TBA and lattice data to differ after some value of r z , at the border of the domain 
where the TBA data are reliable. This is observed in particular for the unsubtracted case, 
where lattice data are affected only by the small uncertainties of the R l — > oo extrapolation. 
In addition, we have observed a large unstability of the lattice data in the region r l large. 
This is very marked for n < in particular, where for R l large enough we have observed 
another pair of extrema in ( |4 . 4| ) as well as related features in the unsubtracted data. An 
example is given for n = — 1 in figure 2. 

5. The boundary case n = — 2 

The case n = —2 is easily solvable. Indeed, the associated Gaussian model at the 
UV fixed point has g = 2, and in the Thirring version, the four-fermion coupling vanishes. 
This is actually the "natural" boundary of the behavior we want to study since the dense 
0(n = —2) model, corresponding formally to g = in ( |3.20| ), is in fact massive. So we 
shall not be able to recover a Gaussian model in the IR, and this example is pathological 
in some respects. 
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That we are dealing with free fermions can also be seen directly from the point of 
view of the 0(n = —2) lattice model. Consider indeed the fermionic integral 

z= j ^# + n e_ ^ + n ex p^^ fe+ ^)]- c 5 * 1 ) 

i (jk) 

By usual rules of Grassman integration this gives rise to an expansion similar to ( |3.1] ) 
with the only difference that "loops" covering a single edge (corresponding to terms like 
tpj~ipk?P'k' l i ; j) are anowe d here. Such a model of loops + dimers j|3| is expected to be in 
the same universality class as (|3.1|) . On the other hand it is exactly solvable by elementary 
techniques since (|5.1| ) is equal to the determinant of the interaction matrix. The critical 
point occurs at (3 C = 1/3 where this matrix restricts to the discrete Laplacian. For (3 close 
to (3 C (mapping between lattice and continuum is straightforward here, so we suppress 
the lattice label) the matrix coincides with A + m 2 with m 2 — 9(/3 c — f3). Therefore the 
massive flow corresponds to a positive m 2 while the massless flow corresponds to negative 
m 2 . These results are in agreement with the value h = h = for the thermal operator of 
the 0(n = —2) model. The sector where non-contractible loops in the spatial direction 
have weight 2 is reproduced by choosing antiperiodic boundary conditions for the fermions 
on the cylinder. In the massive flow the running central charge is immediately obtained 
from [HlHl: 



1 3r 2 / or 1 

1 In- + - + - lE 

7T Z \ r 2 



12 ~ 
+ — / 



7T 

k=l 



y / (2k- l)% 2 + r 2 - (2k - 1)tt - 



2(2k - 1)tt_ 



(5.2) 



where r = mR. In the standard 0(n = —2) model this result is still expected to hold by 
universality, although the correspondence between the mass m and the lattice temperature 
may involve additional logarithmic contributions. 

Consider now the massless direction. The above formula is still expected to hold, 
up to a bulk r 2 (and possible logarithmic corrections) term. The striking fact is then 
the presence of the square root singularities which, beyond the value r 2 = — {2k — l) 2 7r 2 , 
develop an imaginary part. (Similar square-root singularites seem to appear in affine Toda 
theories at negative z |45[]). Recalling the expression for the continuum partition function 
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(including massive states as well as the ground state) of the 0(n = —2) model in that 
particular sector, where T is the length of the cylinder, 

2 + (-) 2 ^ (5.3) 



exp (4tt~) | II 1 + ex P (-!>/(2*-l) 2 T 2 + r 2 ) 



we see for instance beyond r 2 = — n 2 the two largest eigenvalues become equal in modulus. 
Successive branchings occur as r 2 crosses the roots — (2k — 1) 2 tv 2 . It is not clear how to 
define a central charge when the largest eigenvalue becomes complex. The most natural 
candidate is the real part of ( |5.2| ) (the modulus of the largest eigenvalue). Because of 
the square-root term, it is singular with an infinite derivative at each of the points r 2 = 
— (2k — l) 2 7r 2 , and shows wide oscillations moving away from the UV fixed point. This is 
fully confirmed by the numerical study of the lattice 0(n = —2) model ( |3.1| ). 

It is interesting to compare this situation with the square-lattice Ising model. Intro- 
ducing the variable v = tanh/3/2 the locus of zeroes of the partition function is a set of 
two circles depicted in figure 3. The ferromagnetic critical point is the intersection of the 
rightmost circle with the positive real axis. Away from the critical point the Ising parti- 
tion function is a sum of square roots of determinants of the massive Laplacian. The case 
m 2 real and positive corresponds to moving on the real physical axis (in which direction 
depends on the sign of m). The case m 2 real and negative corresponds to moving exactly 
along the circle of singularities. The intersection of the two circles corresponds to m 2 = —4 
which is also the opposite of the maximum eigenvalue of the massless Laplacian on the 
square lattice. 

That the largest eigenvalue of the transfer matrix in a sector can branch and become 
complex for real j3 seems peculiar to n = —2 ||46|| . For any other n G [—2,2] and for R l 
large enough this eigenvalue appears to remain real in our computations. It must be so in 
the r l — > oo limit so that it agrees with the known spectrum in the continuum. 

Although pathological, the case n = —2 is interesting since then the running central 
charge presents an infinite sequence of singularities and oscillations. This is a strong 
indication that for n negative and close enough to n = —2 similar oscillations are present 
(a fact we observed numerically above for n = —1) and maybe also some singularities. 
The region n negative seems more difficult to study because of these features, and in the 
remainder of this paper and the sequel we mainly consider < n < 2. 

As a final comment, we note that the general case should also be amenable by exact 
solution (using the Bethe ansatz) of the Thirring model with imaginary mass. We have 
not done this calculation yet. 
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6. Detailed study of the case n = 0: 
Spontaneous breaking of N=2 supersymmetry 

We now discuss in some detail the case n = 0. It is closely related to the physics of self- 



avoiding walks and also to N=2 supersymmetry [47j. Again, we analyze the ground-state 



energy. It was shown in |§8j that the sine-Gordon model ( |2.5|) is N=2 supersymmetric at 
the coupling where n=0. This will prove a very powerful tool. The fermion number F in 
the N=2 language is equal to half the soliton number in the massive sine-Gordon model. 

The results depend crucially on the boundary conditions around the cylinder |2l[| . 
These are implemented by giving non-contractible loops a weight 2cos7re in the 0(n) 
model, which corresponds to a fugacity for the sine-Gordon solitons of exp±ie7r. In the 
N=2 picture, this corresponds to inserting the operator e l2lTeF . The boundary conditions 
used in sect. 4 have e=0, and in the N=2 language these are called Neveu-Schwarz. These 
give an effective central charge c = 1 in the UV limit. We first discuss the situation with 
twisted boundary conditions, which have a "charge at infinity" e = eo = ±1/2. This gives 
non-contractible loops a weight of zero, and has cjjv = 0. In the N=2 picture this means 
that we insert (— ) F , and is called the Ramond sector. These boundary conditions have 
the advantage that we can derive quantitative information by using the supersymmetry, 
even in the massless phase. 

6.1. The Ramond sector 

a) Preliminary: one-dimensional phase transitions for polymers 

There are several aspects that deserve attention. As a preliminary we recall that for 
lattice cylinders of finite circumference R l , the flow to the dense phase is characterized by 
an infinite sequence of "first-order" phase transitions, taking place at values of f3 > (3 l c 
which depend on R l . These transitions are easily understood by studying the correlation 
functions of L-\eg operators, which are defined by "forcing" a set of L lines to propagate 
on a cylinder (see figure 4) with as usual a weight (3 l per edge. These operators have UV 
conformal dimensions h = h = (9L 2 — 4)/96. The resulting partition function grows like 
the largest eigenvalue A^ ^ of the corresponding transfer matrix. For j3 l < (3 l c self-avoiding 
walks cover a negligible fraction of the lattice, so this eigenvalue must be smaller than 
one. If the total number of possible configurations for this set of L lines of total number of 

(R l ) ( (R l )\~^ 

edges N grows as Q L oc ( fx L ) , they will start occupying a finite fraction of available 

space when j3 l = ^ . In other words, this eigenvalue of the transfer matrix will cross 

ML 
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the previously-dominant trivial eigenvalue, which is always equal to one for any (3 l . The 
polymer density is discontinuous at that point, and this crossing can be described as a 
first-order one-dimensional phase transition H9 . 



To illustrate this point, it is clear that the transfer matrix for propagation of two 
(unoriented) lines on a cylinder of radius R l = 2 is one-by-one, with a single eigenvalue 
equal to 2(f3 1 ) 4 . The latter crosses one for (3 l = 2 -1 / 4 pD[ . 



For fixed L, l//4f ^ ~^ Pi as R l ~ > 00 • Moreover, by the usual finite-size scaling 

"L 

limit, the transitions occur at finite distance from one another. We introduce the notation 



arguments, (l/fi^ ^—(3 l c ) l/ oc 1/R 1 . Therefore in the "blown-up" variables of the continuum 



Jim (l//^ - (3 l c TR l = r l L . (6.1) 
A set of estimates of r l L obtained from lattice computations is given in the following table. 





R l =L 


R l =L+l 


R l =L+2 


R l =L+3 


R l =L+4 




2 


.8101 


.779 


.759 


.746 


.738 — > 


.725 


4 


1.920 


1.788 


1.725 


1.681 


1.653 — ► 


1.6 


6 


3.032 


2.811 


2.700 


2.63 


2.55 — ► 


2.2 


8 


4.147 


3.85 











where the last column is the extrapolation to R l — > 00. They are universal up to an overall 
scale. 

b) Level crossing in the Ramond sector 

By using the quantum-group symmetries discussed in sect. 3, one easily finds that the 
eigenvalues describing L-leg operators are also present in the Ramond sector (e = 1/2) for 
L = 2 mod 4. One can see this explicitly for the value R = 2, where the transfer matrix 
of the 0(n = 0) vertex model in the sector of vanishing spin reads 

{(3 l Y -z(^) 4 

(P 1 ) 4 I . (6.2) 
2(P l ) 2 i 1 / 2 2(P l ) 2 i~ 1 / 2 

The eigenvalues are 1, 2(/3') 4 , 0. Hence at r l 2+4:n a new eigenvalue crosses the value 1, which 
of course is the continuation of the largest eigenvalue in the high-temperature phase. This 
crossing is a first-order phase transition for finite R l and remains so in the R l — ► 00 limit, 
as can easily be checked using scaling or by numerical analysis. 
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Much can be learned about these first-order phase transitions from the N=2 point of 



view. Well-known arguments |5T|] show that the Ramond ground-state energy being zero is 
equivalent to supersymmetry remaining unbroken (for unitary theories) . Since the Witten 
index is not zero here, supersymmetry cannot be broken in the high-temperature phase 
where the model is unitary, and indeed the largest eigenvalue of the transfer matrix is one. 
However, unitarity is lost in the low-temperature phase, and it is possible for a level to cross 
zero energy (the continuation of the Ramond ground state) and become negative. This is 



actually easy to see using the N=2 Landau-Ginzburg action |52| with superpotential 



X 3 + X S + (5(3) 1/2 (X + X), (6.3) 

For (3 > (3 C the potential is not positive (or even real), ensuring that the ground-state 
energy does not have to be zero. We can say that supersymmetry is spontaneously broken 
in the Ramond sector for r l > r 2 , and that at each r l 2+An a new level crosses energy zero. 
Notice however that, since the perturbation does not break supersymmetry explicitly, the 
Witten index remains equal to its UV value, and in fact a pair of levels with opposite 
values of (— ) F cross zero together. 

If one simply measures the ground state energy of the R sector, it is zero up to r\ 
and then non-zero with a discontinuity of the first derivative at that point. However there 
is also a discontinuity of the bulk term, so if one measures running ground-state energies 
after subtracting this bulk term, the final result has a discontinuity at r 2 , as is illustrated 
in figure 5 from lattice data. This is an indication that the analytic continuation done in 
the last section only describes the ground state energy up to r\: since c = in the Ramond 
sector for z > 0, analytic continuation obviously gives c = for all z. Of course all the 
above results apply without reference to the lattice model, and using the variable r. 

This picture can be further elucidated by considering the "index" studied in ]2"2],F23" 



Q = i(R/T)tre- RH F(-) F . i The index Q is obtained by taking the derivative of the 
ground state energy with respect to e at e = 1/2. It has a simple polymer interpretation 



in the massive flow: it is the partition function for a single non-contractible loop |2T| . This 
index has the advantage that we can analytically derive an differential equation for it even 
in the massless phase, as opposed to the ground-state energy in the Neveu-Schwarz sector 
discussed in the last section, where the TBA integral equations can only be numerically 



3 Notice that this Hamiltonian propagates the system in the indirection, as opposed the trans- 
fer matrix we studied above, which is in the T-direction. 
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continued. (If one tries to continue the equations directly, the integrals diverge.) Since all 
of the derivations of |2!| hold for any superpotential, they should hold even in the case 
of a non-real potential like ( |6.3| ) in the dense phase. It was shown in |22|,f23] that in the 
massive phase Q satisfies a differential equation of Painleve III type: we have 

d2u 1 dXJ . 

— + -— = smh[/, (6.4) 
dr z r dr 

with 

r dU 

Near r = mR = the solution we are interested in looks like 

U~^log r -+s + J29nr^ n/3 \ (6.6) 

n 

where exp(s/2) = r(l/3)/(2 2 / 3 r(2/3)) ]53j. As before, we define the variable z = r 4 / 3 . The 
low temperature phase (3 > (3 C corresponds to z negative but real. Defining u = U + in/2, 
(|SlD is then 

d?u 1 du 9 

-r~o H t = ^=cosh it. (6.7) 

dz z z dz 16v— z 



It follows from (|6.6|) that our solution for u and Q is real for negative z. 

For r near r l 2+4n two eigenvalues of the transfer matrix (two levels of the hamiltonian) 
are very close for e = 1/2, and are coupled if e ~ 1/2 by the standard formula of second- 
order perturbation theory. It is easy to show that the numerator is of order e — 1/2 while of 
course the denominator is the difference of energies of the two levels. Therefore we expect 
the derivative with respect to e to exhibit a simple pole at r2+4 n , that is we expect simple 
poles for Q. That Q should become infinite when a level crosses zero energy is clear a priori 
from its definition where no bulk term contribution is subtracted. In terms of polymers 
this pole corresponds to the sudden growth of infinite arms to the non-contractible loop 
(figure 4b). This can be seen explicitly on the lattice for radius R = 2, where for e close 
to 1/2 the eigenvalue that is dominant in the UV reads 

l+(^/3 4 - T ^ i ) 7 r(e-l/2). (6.8) 

Hence Q should exhibit an infinite sequence of simple poles in this massless flow. A 
numerical solution of ( |6.7|) fully confirms this prediction. The locations of the first poles 
(in the variable r) are 
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-^st c^rid grd ^th t^th 

2.95708396 8.61 14.1 19.5 24.8 

Taking ratios to eliminate the lattice-dependent mass scale they agree well with the above 
results of lattice computations. For the reader with this particular interest we give inde- 
pendently in the appendix the results for polymers deduced from this study. 

In fact, we can extract information about these poles analytically. The first pole Z\ 
can be determined from the coefficients g n in ( |6.6| ): as long as there are no other poles 
nearby, the coefficients should be the same as those of the expansion of l/(z — z±) in 
powers of z. The g n can be found analytically by substituting ( |6.6| ) into (|6.4| ); for example, 
gi = — (9/32)2 2 / 3 e _s . Alternatively, one can fit these coefficients to a ratio of polynomials 
as in sect. 4, and find the first pole in this manner (this is why the location of this pole is 
known very accurately). 

It is easy to show that all poles in Q are simple ones. Near a pole, the cosh-Gordon 
equation ( |6.7| ) reduces to the (radial part) of the Liouville equation (i.e. coshu ~ e u /2). 
This can be solved explicitly; the general solution near a pole at z = zq is 

u(z) ~ In — f=t — ttt (6.9) 

Thus all poles in Q oc du/dz are simple. 

We can also see that there are an infinite number of poles on the negative z-axis. We 
show that given one pole, there is another one a finite distance farther along the axis. The 
crucial fact is that in any given region, the right-hand side of (|6.7|) is finite and positive. 
(It may appear more comfortable to suppress the l/\/—z by switching to the variable 
w = (— z) 4 / 3 , where this factor disappears and the right-hand side is > 1.) As we move 
away from the pole in the negative direction, u" and u' are positive. As long as u' / z is 
negative, (|6.7|) requires that u" is positive and finite. Thus u' must decrease at a finite 
rate and eventually hit zero. This is a minimum of u, because on the negative z-axis, u 
has no maxima (only poles): if u' = 0, u" > 0. Thus after the minimum, u must start 
increasing and not turn around. It must increase at a finite rate, so eventually u will get 
large enough so that we can use replace ( |6.7|) with the Liouville equation, and our solution 
is of the form ( |6.9| ). Since u' is negative and must remain negative until we hit a pole, the 
pole in (|6.9|) must fall on the negative axis. Thus we have established that there is another 
pole a finite distance away at larger \z\. 
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In fact, for \z\ large we can show that the distance between a pole at z\ and the next 
pole is ~ Cl^l 1 / 4 , where C is a constant which is about 7. This is done by patching 
together the solution fl6.9Q for u large with the solution for u small (i.e. cosh it ~ 1), which 
is 

c(z) = -z s/2 + Alnz + B. 
4 

As r ^ oo an infinite number of levels crosses the value one. The Witten index 
tr(— 1) F is still equal to 2. However, most levels other than the ones of vanishing energy 
have a bulk term that becomes infinite in that limit. If we subtract this bulk term, we 
have to set the Witten index to zero for consistency. This is in agreement with the study 
of the IR fixed point described by an ?y£ system. We have here two theories related by 
a N=2 preserving flow which have different Witten indices, so they must be considered 
"at infinite distance" from each other p4||. This sheds some doubt about the possibility 



of connecting them smoothly by a TBA-like description, an issue which we will discuss 
further in the next paper. 

Besides the crossing of eigenvalues at r l 2+An just explained by quantum-group sym- 
metries, we have observed additional crossings at values that tend toward r l A+An for large 
lattice systems. The matrix elements are such that these crossings do not produce a pole 
in Q. 

6.2. The Neveu-Schwarz sector 
a) Observations 

This is the sector discussed in sect. 4, where solitons for the massive flow in the sine- 
Gordon model have fugacity one, and non-contractible loops in the 0(n) vertex model are 
given a weight 2. Unfortunately, there is no index in this sector and the TBA continuation 
of sect. 4 breaks down after the first maximum. We here discuss only the lattice results. As 
in the large-t expansion, c = 1 is observed both in the UV and IR. A remarkable feature, 
very different from the massive study, is the very slow stabilization of the curve's shape as 
the lattice size is increased, and one may wonder if this hides further aspects, unobservable 
at the sizes we can study. However, a clear result which generalizes to all values of n > —2 
is that the largest eigenvalue of the transfer matrix remains real and non-degenerate at 
all values of r . For the moment we start by doing some "numerology". Comparison of 

, ,2 

estimates for r L with figure 6 shows that peaks of the derivative d *j~ r i c y occur roughly at 
r\ and r\. These coincidences are confirmed by a more detailed study, and they hold 
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within our error bars. At r l = r 2 there occured a phase transition in the Ramond sector, 
associated with a discontinuity of the running effective central charge. It is not surprising 
that r\ emerges as well in the study of the NS sector. Even though we are a finite distance 
from the UV fixed point in the blown- up variable r , we are still infinitesimally close to it, 
and the polymers present in the ground state of this sector still occupy an infinitesimally 
small fraction of the available space. This does not change the fact that at r l 2 two-leg 
diagrams become infinite. One may actually wonder whether r\ (as well as r l 2n ) could be 
a point of phase transition in the sine-Gordon model itself, independent of the sector; we 
shall discuss this further. 

We have also observed that the two (first) extrema of the curve c(r ) occur roughly 
at the values r\ and r l 3 . Study of the present data indicates that they do not exactly take 
place there, although as always this may be a finite-size effect. We do not expect any 
special feature to occur in the Ramond or Neveu-Schwarz sector at values r l 1+2n . This 
is because the transfer matrix we are considering acts on spin zero in the vertex model 
language, so only even numbers of lines can propagate through the system. 

b) Possible singularities 

The evidence seems to exclude a discontinuity in c(mR). As far as its derivatives are 



concerned, they are difficult to study numerically. The lattice-subtracted quantity (|4.4|) is 
not very reliable for determining derivatives, nor are any of the procedures we have used 
to subtract the bulk term. We can of course study derivatives of unsubtracted quantities, 
but then their convergence as L increases is affected by c and the bulk term, and this 
mixture is difficult to disentangle. The safest quantity to study seems to be the second 
derivative with respect to r of the full ground-state energy, including the bulk contribution. 
In the large R l limit, the bulk term should simply contribute an additive constant. Results 
are shown in figure 7. There is a spectacular peak at values that converge to r 2 . The 
convergence past this value is however quite slow, and makes the quantitative analysis of 
this peak difficult. We have done the same analysis for the Ramond sector, where we know 
that we should observe a delta function. The numerical results look virtually identical 
to these past r 2 . We therefore suggest that there may be a discontinuity of the first 
derivative of c in the Neveu-Schwatz sector at r 2 . The analysis of shows similar but 
less marked properties. We could not get any information about other r 2n . However if 
there is indeed a singularity for r 2 and it seems natural there is one at all r 2n . Recall 
that in the polymer language, at r = r<i n it becomes possible for 2n lines to propagate 
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along the cylinder in time direction, and that this interferes with the characteristic feature 
of the ground state in the UV limit, the domination by non-contractible loops in space 
direction (with weight 2). 

To investigate further the possibility of a singularity at r\ we studied in detail the 
structure of the transfer matrix in the Neveu-Schwarz sector for complex (3 . For finite 
lattice size R l , the largest eigenvalue does not cross for (3 real. One can however locate 
square root branch points in the complex plane, a pair of which moves (in the blown-up 
variable r l ) towards r\. The situation here is somewhat simpler than for ordinary critical 
phenomena, since in the variable r l the infinite number of branch points that converge to 
the critical inverse temperature (3 l c are now at finite distance from each other. Hence we 
have to consider the limit point of a single pair of branch points, which is expected to 
exhibit a singularity of the form x k \x\ (see [55] for more detailed discussion). This agrees 



with the above expectation of k = 0. However a singularity of the largest eigenvalue is not 
enough here; we need to make sure that it occurs in its scaling part. We unfortunately did 
not get definitive numerical evidence for this. 

6.3. Miscellaneous 

We have two additional observations. At the UV and IR fixed points, we have h = 
gm 2 /4 in the NS sector. Thus we can define the running coupling constant g{r l ) by 
measuring the gap between the ground state of the NS sector and the ground state of the 
sector where two non-contractible loops propagate (m=2). The result for available sizes is 
a monotonic curve with a derivative maximum at a value close to g = 1, close also to r\ 
(these coincidences do not seem to be exact). 

We have also investigated whether the model passes "close" to an intermediate fixed 
point by studying the spectrum of excited states and looking for the tower structure char- 
acteristic of a conformal field theory. No such structure has been observed. 



7. Conclusions 

We have explained why n negative was expected to exhibit even more complicated 
properties, and we shall not discuss it further. For n > 0, the key features of the large-t 
expansion are reproduced by the various numerical studies. In particular, we emphasize 
that despite the complex nature of the interaction, the ground state energy of the sine- 
Gordon model (except for n = —2) remains real all along the flow. We have also checked 
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that the levels of the sine-Gordon model that contribute to the minimal partition coincide 
with predictions from the massless TBA for minimal models. We have presented a detailed 
study of the case n = because supersymmetry allowed us to do so. This case is certainly 
a little marginal, since for instance the associated c = "minimal model" has no degrees 
of freedom and a partition function equal to one, and so is clearly "infinitely far" from the 
non-minimal one whose partition function grows exponentially. This case presents strong 
evidence for genuine singularities in the flow whose physical interpretation is rather clear 
from polymer point of view. The natural questions to ask now are 

1) are there really singularities in the flow for n = 0? 

2) if there are indeed singularities for n = 0, is it a marginal case or do all < n 
values exhibit the same behavior? 

Certainly for < n < 2 the situation is somewhat different than for n = since the 
free energies of minimal and sine-Gordon models are then the same (level crossings occur 
in some sectors, but only in finite numbers). 

3) is there an exact scattering theory for massless flow within sine-Gordon and between 
successive minimal models, at least for some values of < n? 

The study of these questions and partial answers are presented in the next paper. 
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Appendix A. Polymer results 

We summarize here the polymer results obtained in this paper for the reader with this 
specific interest only. Consider the total number of configurations for a set of L polymer 
lines of total number of edges M crossing the cylinder of radius R l (which we denote simply 
by R here). One has 

rt L R) oc [^ ] ) N , M»l. (A.l) 

The connectivity constant //^ is such that for every L, fj,^ — > [i as R — ► oo where \i is 
the usual lattice-dependent "bulk" connectivity constant [|56) . Using finite-size scaling one 
has /i^ — /i oc J2 _1 / I/ , where v = 3/4. Define the numbers (where j3 c = 

linn^l//^ - PcY& = r l L . (A.2) 
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These numbers are universal up to a lattice-dependent overall scale (that is r^jr^i is 
universal). We have found that the particular numbers r2+4 n are in the same ratio as the 
(simple) rfi 1 poles of the function Q given by ( |6.5|) ; Q is the derivative of a solution of the 
Painleve III equation fl6.4j ), with boundary condition (|6.6|) . The locations of these poles 
do not have simple expressions. The first few numerical values are 

-^st ^nd grd ^th r^th 

2.95708396 8.61 14.1 19.5 24.8 

We have also found that the distance r 2 +4n — r_ 2 +4 n behaves as Cr 4n for n large, where 
C is a constant close to C ~ 5. 

That Q should have a singularity can also be established on physical grounds using 



the polymer model. Indeed recall that in [21] we showed that 



oo 

4/3 



Q{r) = / e" (r) x f(x)dx. (A.3) 



o 

where the number of configurations of a single non-contractible loop of length M on the 
cylinder reads in the scaling limit 

u R = TR-y*^sf(^- 3 ). (A.4) 

and we have set 

while s is some lattice-dependent rescaling factor. The dominant behavior of the function 
/ follows from the loss of entropy due to confinement that is f(y) ~ y a exp(— cy) where 
a, c are constants to be determined. Plugging this result back into the Laplace transform 
( |A.3| ) we see immediately that c = (r2) 4//3 and a = 0. That the connectivity constant for 
a non-contractible loop and for two propagating legs (L = 2) should be the same appears 
reasonable, and it is proven by the quantum-group arguments in the text. 

Some of these results could likely be extended to the other number of legs, in particular 
odd ones. However the connection with the PHI equation, which has been our main 
analytical tool here, is possible only for a number of legs equal to two modulo four, which 
is characteristic of the Ramond sector in the N=2 point of view. 
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Figure Captions 



Figure 1: Data for the running central charge c(mR) for n = V2 obtained from the 
lattice model and analytic continuation of massive TBA 

Figure 2: Data for the running central charge at n = — 1 obtained from the lattice 
model. Observe the appearance of additional extrema in the IR as R l increases. 

Figure 3: The locus of zeroes for the Ising model partition function. Perturbation by 
a purely imaginary mass moves along the circles. 

Figure 4a: A set of L polymer lines propagating on the cylinder 

Figure 4b: Growth of an infinite arm at r 2 

Figure 5: Discontinuity of the running subtracted central charge in the Ramond sector 
(n = 0) at r 2 

Figure 6: Running central charge c{mR) for n = obtained from the lattice 
Figure 7: Second derivative of c(r = mR) with respect to r for n = 0. 
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Fig. 1. Comparison of TBA and lattice data for n=/2 




Fig. 2. c(r) for n=- 1 




Fig. 4a. Two polymers propagating together 




Fig. 4b. Growth of an infinite arm 




Fig. 5. c eff in Ramond sector for n=0 



